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ABSTRACT 

We use numerical hydrodynamical simulations of early structure formation in a ACDM universe to 
investigate the spin temperature and 21cm brightness of the diffuse intergalactic medium (IGM) prior 
to the epoch of cosmic reionization, at z < 20. In the absence of any radiative heating, collisions 
between neutral hydrogen atoms can efficiently decouple the spin temperature from the CMB only in 
dense minihalos and filaments. Shock heated gas shines in emission, surrounded by cooler gas visible 
in absorption. In the case of a warm, mostly neutral IGM, produced here by X-ray emission from an 
early miniquasar, the 21cm signal is strongly enhanced. Even slightly overdense filaments now shine 
in emission against the CMB, possibly allowing future radio arrays to probe the distribution of neutral 
hydrogen before reionization. 

Subject headings: cosmology: theory - diffuse radiation - galaxies: evolution - intergalactic medium 



1. INTRODUCTION 

It has long been known that neutral hydrogen in the 
diffuse intergalactic medium (IGM) and in gravitation- 
ally collapsed structures may be directly detectable at 
frequencies corresponding to the redshifted 21cm line 
(Field 1959; Sunyaev & Zel'dovich 1975; Hogan k Rees 
1979). The emission or absorption of 21cm photons from 
neutral gas is governed by the spin temperature T5, de- 
fined as ni/uQ = 3 exp(— T^/Tg). Here uq and ni are 
the number densities of atoms in the singlet and triplet 
n = 1 hyperfine levels, and = 0.068 K is the tem- 
perature corresponding to the energy difference between 
the levels. To produce an absorption or emission signa- 
ture against the cosmic microwave background (CMB), 
the spin temperature must differ from the temperature 
of the CMB, TcMB = 2.73 (1 z) K. At 30 < z < 200, 
prior to the appearance of nonlinear baryonic objects, 
the IGM cools adiabatically faster than the CMB, spin- 
exchange collisions between hydrogen atoms couple Ts to 
the kinetic temperature Tk of the cold gas, and cosmic 
hydrogen can be observed in absorption (Scott & Rees 
1990; Loeb & Zaldarriaga 2004). At lower redshifts, the 
Hubble expansion rarefies the gas and makes collisions 
inefficient: the spin states go into equilibrium with the 
radiation, and as Ts approaches Tcmb the 21cm signal 
diminishes. It is the first luminous sources that make 
uncoUapsed gas in the universe shine again in 21cm, by 
mixing the spin states either via Lya scattering or via 
enhanced free electron-atom collisions (e.g. Madau et al. 
1997, hereafter MMR; Tozzi et al. 2000; Ciardi & Madau 
2003; Gnedin & Shaver 2004; Nusser 2005). 

While the atomic physics of the 21cm transition is well 
understood in the cosmological context, exact calcula- 
tions of the radio signal expected during the era between 
the collapse of the first baryonic structures and the epoch 
of complete reionization have been difficult to obtain. 
This is because the "21cm radiation efficiency" depends 
on spin temperature, gas overdensity, hydrogen neutral 
fraction, and line-of-sight peculiar velocity (e.g. MMR; 

Electronic address: |mqk@ucolick.org| 



Bharadwaj & Ah 2004). When Ts = Tk, the visibil- 
ity of the IGM at 21cm revolves around the quantity 
{Tk - Tcmb)/Tk. If Tk < Tcmb, the IGM wiU appear 
in absorption against the CMB; in the opposite case it 
will appear in emission. To determine the kinetic tem- 
perature of the IGM during the formation of the first 
sources, one needs a careful treatment of the relevant 
heating mechanisms such as photoionization and shock 
heating. In addition to the signal produced by the "cos- 
mic web" , minihalos with virial temperatures of a few 
thousand kelvins form in abundance at high redshift, and 
arc sufficiently hot and dense to emit collisionally-excitcd 
21 cm radiation (Iliev et al. 2002). 

In this Letter, we use fully 3D Eulerian cosmological 
hydrodynamical simulations of early structure formation 
to make detailed predictions of the thermal history, spin 
temperature, and 21cm brightness of neutral hydrogen 
in the pre-reionization era at z < 20. The adaptive mesh 
refinement (AMR) technique allows us to resolve, with 
progressively finer resolution, the radio signal expected 
from the densest parts of the IGM and from miniha- 
los, and study the impact on 21cm spectral features of 
sources of X-ray radiation that may start shining before 
a univeral Lya and Lyman-continuum background is ac- 
tually established. 

2. SPIN-KINETIC TEMPERATURE COUPLING 

In the quasi-static approximation for the populations 
of the hyperfine levels, and in the absence of radio 
sources, the H I spin temperature at a given redshift is 
a weighted mean between Tk and Tcmb, 

T"* + Tcmb + IjTk 



1 + y 

The coupling efficiency y is the sum of three terms, 
T, 

y 



AT, 



(Ch +Ce+ Cp), 



K 

where A = 2.85 x 10^^^ s^^ is the spontaneous emis- 
sion rate and Ch, Ce, and Cp are the de-excitation rates 
of the triplet due to collisions with neutral atoms, elec- 
trons, and protons. A fourth term must be added in the 



(1) 



(2) 



2 



Spin temperature and 21cm brightness of IGM 



presence of ambient Lya radiation, as intermediate tran- 
sitions to the 2p level can mix the spin states and couple 
Ts to Tk, the "Wouthuysen-Field" effect (Wouthuysen 
1952; Field 1958; Hirata 2005). 

The H-H collision term can be written as Ch = tly^hl, 
where k is the effective single-atom rate coefficient re- 
cently tabulated by Zygelman (2005) for temperatures 
in the range I < Tk < 300 K. For 300 < Tk < 1000 K 
we have used the values tabulated by Allison & Dalgarno 
(1969) multiplied by a factor of 4/3 (as recommended by 
Zygelman 2005). These rates can be fitted by a sim- 
ple function, k = 3.1 x 10"" T^-^" exp(-32/Tif) cm^ 
s~-^, which is very accurate (to < 0.5%) in the range 
10 < Tr: < 10^ K. This fit also extrapolates well to the 
Tk = 3000 K and Tk = lO^'K coefficients (x4/3) listed 
in Field (1958). For the e-H collision term, Ce = fi-eTe, we 
have used the functional fit of Liszt (2001) to the down- 
ward rates of Smith (1966), log(7e/cm3 s^^) = -9.607 
0.51og(Tif)exp[-(logTif)'' Vl800] for Tk < lO^K, and 
-fe{TK > 10'* K) = 7e(10''K). The rate coefficient for 
proton de-excitation is just 3.2 times larger than that 
for neutral atoms at Tk > 30K (Smith 1966). We have 
used 7p = 3.2 k (with Cp = ripjp): excitation by pro- 
tons is typically unimportant as it is much weaker than 
that by electrons at the same temperature. In the ab- 
sence of Lya pumping, to unlock the spin temperature 
of a neutral medium with Tk = 500 K from the CMB 
requires a collision rate Ch > ATqmb/T^, or a baryon 
overdensity Sb > 5 [(1 + z)/20]^^. Not only dense gas 
within virialized minihalos but also intergalactic filamen- 
tary structure heated by adiabatic compression or shock 
heating may then be observable in 21cm emission. We 
shall see in the next sections how a population of "mini- 
quasars" turning on at early stages (e.g. Madau et al. 
2004; Ricotti et al. 2005) will make even the low-density 
IGM visible in 21cm emission as structure develops in 
the pre-reionization era. 

3. NUMERICAL SIMULATIONS 

The simulations have been performed using a modified 
version of enzo, an AMR, grid-based hybrid (hydro+N- 
body) cosmological code developed by Bryan & Norman 
(sec http://cosmos.ucsd.edu/cnzo/). Details of the sim- 
ulation setup are discussed in Kulilen & Madau (2005, 
hereafter KM05). We follow the non-equilibrium, nine- 
species (H, H+, H-, e", He, He+, He++, H2, and H2+) 
chemistry of primordial self-gravitating gas, including ra- 
diative losses from atomic and molecular line cooling, sec- 
ondary ionizations and Compton heating by X-rays, and 
Compton cooling by the CMB. These simulations dif- 
fer from those described in KM05 only in the use of the 
H2 cooling function of Galli & Palla (1998) (instead of 
Lepp & ShuU 1983) and of a new solver that couples the 
nine-species rate equations with the radiative gas cooling 
equations. The AMR is restricted to the inner 0.5 Mpc 
of a box of i = 1 Mpc (comoving) on a side: this region 
is allowed to dynamically refine to a total of 8 levels on 
a 128'^ top grid, resulting in a maximum resolution of 
30 pc (comoving).^ While the L = 1 Mpc box was se- 
lected to be a representative patch of the universe (with 

^ The DM density field is also sampled with 128^ particles, lead- 
ing to a mass resolution of moM = 2000 Mq. This ensures that ha- 
los above the cosmological Jeans mass are well resolved at ^ < 20. 



mean matter overdensity S ~ p/p ~ 1), the inner 0.5 
Mpc region is characterized by S = 1.37, corresponding 
to a 1.2 — cr fiuctuation. 

The two-dimensional distribution of gas overdensity 
and spin temperature at z = 17.5 is shown in Figure ^ 
for a simulation with no radiation sources ( "NoBH" ) and 
one ("PL") in which a miniquasar powered by a 150 Mq 
black hole turns on at redshift 21 within a host halo of 
mass 2 x 10^ M©. The miniquasar shines at the Edding- 
ton rate and emits X-ray radiation with a power-law en- 
ergy spectrum oc in the range 0.2—10 keV. The black 
hole mass grows exponentially to 400 Mq by z = 17.5. 
The color coding in this phase diagram indicates the frac- 
tion of the simulated volume at a given (Sb, Ts). In both 
runs we have assumed no Lya pumping, so that the vis- 
ibility of hydrogen at 21cm is entirely determined by col- 
lisions. Only gas with neutral fraction > 90% is shown 
in the figure. The low-density IGM in the NoBH run 
lies on the yellow Ts = Tcmb = 50.4 K line: this is gas 
cooled by the Hubble expansion to Tk ^ Tcmb that 
cannot unlock its spin states from the CMB and there- 
fore remains invisible. At overdensities between a few 
and ^ 200, H-H collisions become efficient and adiabatic 
compression and shocks from structure formation heats 
up the medium well above the radiation temperature. 
The coupling coefficient at this epoch is y ^ SbT^^'^"^: 
gas in this regime has Tcmb < Ts ^ vTk < Tk and 
appears in emission against the CMB (red and green 
swath). Some residual hydrogen with overdensity up to 
a few tens, however, is still colder than the CMB, and is 
detectable in absorption. At higher densities, y 3> 1 and 
Ts — * Tk- the blue cooling branch follows the evolution- 
ary tracks in the kinetic temperature-density plane of gas 
shock heated to virial values, Tk = 2000—10'* K, which is 
subsequently cooling down to ^ 100 K because of H2 line 
emission. The volume and mass fractions of gas with 
Ts > Tcmb within the simulation box are fv ~ 0.002 
and /m = 0.068, respectively. The latter is comparable 
to the amount of gas in the shocked phase estimated by 
Furlanetto & Loeb (2004, their fig. 8) using an extension 
of the Press-Schechter formalism. 




Fig. 1. — Two-dimensional distribution of spin temperature ver- 
sus baryonic overdensity at 2 = 17.5. The color coding indicates 
the fraction of the simulated volume at a given {&\, , Tg ) . Le/i; 
NoBH run. The volume and mass-averaged spin temperatures are 
48.6 and 67.5 K, respectively. Right: PL run. Only gas with 
neutral fraction > 90% is shown in the figure. The volume and 
mass-averaged spin temperatures arc 82.6 and 138.0 K, respec- 
tively. The red horizontal line marks the temperature of the CMB 
at that redshift. 

The effect of the miniquasar on the spin temperature 
is clearly seen in the right panel of Figure ^ X-ray 



Kuhlen, Madau. & Montgomery 



3 



radiation drives the volume-averaged temperature and 
electron fraction (xe) within the simulation box from 
(8K,1.4 X 10-4) to (2800K,0.03), therefore producing 
a warm, weakly ionized medium (KM05). The H-H col- 
lision term for spin exchange in the low-density IGM in- 
creases on the average by a factor 350° '^^ ~ 8, while 
the e-H collision term grows to Ce ^ 0.5Ch- Gas with 
{Si,,TK,Xe) = (1, 2800 K, 0.03) has coupHng efficiency 
y = 0.008, spin temperature Tg = 73 K > Tcmb, and can 
now be detected in emission against the CMB. Within 
150 comoving kpc from the source, the volume-averaged 
electron fraction rises above 10%, and e-H collisions dom- 
inates the coupling. 

4. RADIO SIGNAL 

A beam of 21cm radiation passing through a neutral 
hydrogen patch having optical depth r and spin temper- 
ature Ts causes absorption and induces emission. In the 
comoving frame of the patch, the radiative transfer equa- 
tion yields for the brightness temperature through the 
region: Tb = Tqubs''^ + Ts{l — e^"^). We have used our 
numerical simulations to perform 21cm radiation trans- 
port calculations, including the effect of peculiar veloci- 
ties and local changes in spin temperature, gas density, 
and neutral hydrogen fraction, as follows. The equation 
of radiative transfer has been discrctized according to 



Tb,i^o — Tcmb e 



(3) 



where Tg is the spin temperature of the i cell and 
is the line-center 21cm frequency as measured in the box 
plane where the miniquasar is located. The i-sums are 
over all cells along an orthogonal line-of-sight through 
the box, while the j-sum is over all cells in front of the 
z**^ cell. The optical depth of a cell with H I column iV^j 
is then given by 



At,! 



(4) 



where Aiy — v^Vi/c to account for the bulk motion as- 
sociated with peculiar and Hubble flow velocities, and 
0(Az^) is the normalized line profile. The total line-of- 
sight bulk velocity of the i^^ cell relative to the mini- 
quasar plane at proper distance can be written as 
V, = v^^p+H{z)n, where H{z) = Hq [nM{l+zy^+^AY^^ ■ 
In our simulations the typical line-of-sight component 
of the peculiar motion dominates over the Hubble ex- 
pansion on scales below 40 comoving kpc. The line 
profile is entirely determined by thermal broadening, 
so that (f){Au) ~ (Ai^uy^)-^ exp[— (Ai//Az/£))^], where 
Avj) = vq^^ 2kTK / (mpC^) is the Doppler width. The 
final step in the process is to propagate the 21cm radi- 
ation from the simulation box at redshift z = 17.5 to 
the Earth. Neglecting foreground emission and absorp- 
tion, the observed flux at frequency i'o/(l + z) can be 
simply expressed by the differential brightness tempera- 
ture against the CMB as ST^ = (Tb - Tcmb)(1 + ^)"^- 
We have calculated 5Tb maps consisting of 1024^ pixels, 
corresponding to an angular resolution of ~ 0.01" and a 
frequency resolution of 21 Hz. The latter is three orders 
of magnitude smaller than the Hubble flow across the 
entire box. The resulting 21cm radio signal is shown in 
Figure 12 (lower panels), together with an image of the 



projeciec! hydrogen spin temperature (upper panels): the 
latter highlights the abundance of structure within our 
simulation box on scales up to hundreds of kpc. Due to 
Hubble and peculiar velocity shifts not all of this struc- 
ture contributes to the STb map. In the NoBH simula- 
tion, coherent features in the IGM can be discerned in 
emission (Ts > Icmb, 5Tb > 0): this filamentary shock- 
heated structure is typically surrounded by mildly over- 
dense gas that is still colder than the CMB and appears 
in absorption {Ts < Tcmb, 5Tb < 0). The covering factor 
of material with STb < — 10 mK is 1.7%, comparable to 
that of material with STb > +10 mK: only about 1% of 
the pixels are brighter than -1-40 mK. While low-density 
gas (black color in the left-lower panel) remains invisible 
against the CMB in the NoBH run, the entire box glows 
in 21cm after being irradiated by X-rays. The fraction of 
sky emitting with STb > {+10, +20, +30, +50, +100) mK 
is now (0.57,0.31,0.19,0.1,0.035). 

5. CONCLUSIONS 

We have performed high-resolution AMR cosmologi- 
cal hydrodynamic simulations of early structure forma- 
tion in a ACDM universe to investigate the spin tem- 
perature and 21cm brightness of the IGM prior to the 
epoch of cosmic reionization. In the adopted cosmology, 
the simulated box size corresponds to an angular scale 
A0 = 0.32'(L/Mpc) at z = 17.5. We have not tried to 
match this scale to the resolution/sensitivity of planned 
low-frequency arrays such as LOFAR, MWA, PAST, and 
SKA. Rather, our main goal here was to locate and re- 
solve some of the physics that determines the observ- 
ability of 21cm spectral features in the pre-reionization 
era. We have shown that, even in the absence of ex- 
ternal heating sources, spin exchange by H-H collisions 
can make filamentary structures in the IGM (heated by 
adiabatic compression or shock heating) observable in 
21cm emission at redshifts z < 20. Some cold gas with 
overdensities in the range 5-100 is still detectable in ab- 
sorption at a level of STb — lOmK, with a signal that 
grows as Tg^ and covers a few percent of the sky. X- 
ray radiation from miniquasars preheats the IGM to a 
few thousand kelvins and increases the electron fraction: 
this boosts both the H-H and the e-H coUisional cou- 
pling between Ts and Tk, making even low-density gas 
visible in 21cm emission well before the universe is sig- 
nificantly reionized. Any absorption signal has disap- 
peared, and as much as 30% of the sky is now shining 
with 5Tb > +20 mK. As pointed out by Nusser (2005), 
the enhanced e-H coupling makes the spin temperature 
very sensitive to the free-electron fraction: the latter is 
also a tracer of the H2 molecular fraction in the IGM. 
Note that when Ts 2> Tcmb , the radio signal saturates in 
emission. 21cm saturation can readily be inferred for the 
dark-blue and reddish portions of the upper-right image 
in Figure El The brightness temperature through these 
regions would then remain approximately the same even 
in the presence of background Lya radiation above the 
"thermalization" value for efficient hyperfine level mixing 
(corresponding at these epochs to ~ 0.1 Lya photons per 
hydrogen atom, MMR). Of course, in this case even the 
low-density IGM would be lit up at 21cm, and the overall 
contrast in the radio signal would actually decrease. 

Our NoBH 5Tb maps seem to be in qualitative agree- 
ment with the results of Ahn et al. (2005), exhibiting 
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Fig. 2. — Projected (mass-weighted) spin temperature (upper panels, logarithmic scale) and 21cm differential brightness temperature 
(lower panels, composite linear scale) in a 0.5 Mpc simulation box for runs "NoBH" (left) and "PL" (right) at z = 17.5. The location of 
the miniquasar is indicated by crosses in the right panels. 
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both emission from minihalos and dense filaments and 
absorption from the cold IGM. Gas belonging to virial- 
ized minihalos (i.e. having baryonic overdensity Sb > 70, 
see KM05) contributes about 60% of the total emission 
signal while most of the absorbing signal (light blue color 
in the left-lower panel of Fig. 13) is due to gas with (5b < 10 
still left below Tcmb by adiabatic heating. In the PL case 
instead, the signal is entirely in emission and is domi- 
nated by the warm, mostly neutral, diffuse IGM. Finally, 
we point out that, on the simulated scales, we find lit- 
tle evidence for a systematic enhancement of the radio 
signal due to velocity compression from baryonic infall 
(Barkana & Loeb 2005). Regions that are able to coUi- 
sionally decouple their spin temperature from the CMB 



are non-linear density perturbations, for which redshift- 
space smearing ( "fingers of god" ) outweighs the enhance- 
ment from the Kaiser effect. Peculiar velocities do, how- 
ever, Doppler shift features into and out of the STt maps. 
We are running a suite of hydrodynamic simulations in 
larger cosmological volumes, which will allow predictions 
at the spatial and spectral resolution of future planned 
radio arrays. 
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